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0\ : Abstract 



The dynamics of an interface between the normal and superconducting phases 

o ■ 

CN ' under nonstationary external conditions is studied within the framework of 

the time-dependent Ginzburg-Landau equations of superconductivity, modi- 
' fied to include thermal fluctuations. An equation of motion for the interface 

o ■ 

' is derived in two steps. First, the method of matched asymptotic expansions 

m 



is used to derive a diffusion equation for the magnetic field in the normal 



' phase, with nonlinear boundary conditions at the interface. These boundary 

a" 

conditions are a continuity equation which relates the gradient of the field at 



the interface to the normal velocity of the interface, and a modified Gibbs- 
Thomson boundary condition for the field at the interface. Second, the bound- 



^ ' ary integral method is used to integrate out the magnetic field in favor of an 



equation of motion for the interface. This equation of motion, which is highly 
nonlinear and nonlocal, exhibits a diffusive instability (the Mullins-Sekerka 
instability) when the superconducting phase expands into the normal phase 
(i.e., when the external field is reduced below the critical field). In the limit of 
infinite diffusion constant the equation of motion becomes local in time, and 
can be derived variationally from a static energy functional which includes 
the bulk free energy difference between the two phases, the interfacial energy, 
and a long range self-interaction of the interface of the Biot-Savart form. In 
this limit the dynamics is identical to the interfacial dynamics of ferrofluid 
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domains recently proposed by S. A. Langer et al. [Phys. Rev. A 46, 4894 
(1992)]. As shown by these authors, the Biot-Savart interaction leads to me- 
chanical instabilities of the interface, resulting in highly branched labyrinthine 
patterns. The application of these ideas to the study of labyrinthine patterns 

in the intermediate state of type-I superconductors is briefly discussed. 
74.55.+h, 74.60.-w, 05.70.Ln 
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Typeset using REVT^ 



I. INTRODUCTION 



Over the last decade there has evolved a fairly complete understanding of the physics of 
several types of nonequilibrium growth patterns, such as the dendritic growth of solidifying 
systems (e.g., "snowflakes" ) or the fingered growth which occurs at the interface of driven 
immiscible fluids (e.g., viscous fingers in Hele-Shaw cells); for reviews of the theoretical 
situation, see Kessler et al. [jT| and Langer 0. The similar patterns which grow in these 
ostensibly different physical systems are the consequence of a competition between a dynamic 
instability (the MuUins-Sekerka instability [^] for dendritic growth) which promotes the 
growth of a highly ramified interface, and surface tension, which favors a smooth interface. 
One of the important theoretical insights which has emerged from this work is that surface 
tension anisotropy plays a crucial role in determining the morphology of the pattern 

It has recently been shown that the process of magnetic flux expulsion in type-I su- 
perconductors subjected to a magnetic field quench share many features with these other 
pattern forming systems 0,|5[. In particular, as a superconducting nucleus grows, the ex- 
pelled flux generates eddy currents in the normal phase; the magnetic field h in the normal 
phase therefore satisfies a diffusion equation, 

dth = DV^h, (1.1) 

where the magnetic diffusion constant D = l/Ana, with a the normal state conductivity 
(we will set c = 1). By applying Maxwell's equations to the interface itself, and noting that 
the electric and magnetic fields both vanish in the superconducting phase, we arrive at a 
continuity equation for the field at the interface 

{V xh) xii\i = -Dvnhi, h -111^ = 0, (1.2) 

where n is the unit normal at the interface, directed toward the normal phase, and f„ is 
the normal velocity at the interface. Finally, the field at the interface should equal the 
superconducting critical field He, with curvature corrections: 
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ihi. = |i - 1^ wud) + <sm /c} , (1.3) 

where o"ns(6') is the surface tension of the normal/superconducting interface (not to be con- 
fused with the conductivity), which depends on the angle 6 with respect to the crystal axes (it 
would depend on two angles 61^2 in three dimensions), and K. is the curvature of the interface 
(or the sum of the principle curvatures in three dimensions) • Similar equations (with- 

out the surface tension) were used by Pippard P] and Lifshitz to discuss the growth of the 
normal phase into the superconducting phase, which is a dynamically stable process. The 
analogy with the solidification problem is apparent when the magnetic field is identified with 
the temperature of the liquid, and the magnetic diffusion constant with the thermal diffusiv- 
ity; the continuity condition, Eq. (|1.2| ), is replaced with the continuity condition for the heat 
generated at the solidifying interface, and Eq. ( |1.3|) is a modified form of the Gibbs-Thomson 
boundary condition Due to this formal similarity between the equations describing 

solidification and those describing the kinetics of the normal/superconducting transition, it 
was predicted that the growth of a superconducting nucleus into the "supercooled" normal 
state would be dynamically unstable Such instabilities were observed in numerical 

studies of the time-dependent Ginzburg-Landau (TDGL) equations of superconductivity for 
propagating interfaces; see e.g.. Fig. 1 of Ref. @]. 

In this paper we will connect the TDGL equations to the "sharp-interface" equations 
discussed above. In Sec. II the TDGL equations, including fluctuations, are presented. In 
Sec. Ill we use the method of matched asymptotic expansions to reduce these equations to a 
set of sharp interface equations, generalized to include fluctuations. In Sec. IV we formally 
reduce the sharp interface equations to an equation of motion for the interface itself, and 
discuss some of the features of this equation of motion. In particular, we will show that in the 
limit that the diffusion constant D ^ 00, the interface equation of motion becomes local, 
and can be derived from an interfacial energy functional. This functional contains three 
terms: the bulk free energy difference between the normal and superconducting phases, the 
interfacial surface energy, and a Biot-Savart interaction of the interface with itself. When 
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written in this form, the equation of motion is identical to a phenomenological equation of 
motion for the dynamics of two dimensional domains of ferrofiuids, recently proposed by 
Langer et al. |Ty]. Appendix A reviews some properties of the surface tension, along with 
a new result on the behavior of the surface tension near the critical value oi k = 
Appendix B includes some details on the calculation of the kinetic coefficient which appears 
in the equation of motion. In Appendix C the noise correlations are calculated and shown to 
satisfy the fluctuation-dissipation theorem. Some useful definitions and results concerning 
the differential geometry of curves in two dimensions are provided in Appendix D. 



II. THE TIME-DEPENDENT GINZBURG-LANDAU MODEL 

The time-dependent Ginzburg-Landau (TDGL) model consists of equations of motion 
for the complex superconducting order parameter ip, the magnetic vector potential a, and 
the scalar potential cf) [|rT|. The origin and validity of the equations have been extensively 
discussed elsewhere ||T2|,|13|. In conventional units, these equations are 



h^{dt + ij(t>)i^ = |^(V - ij^fij + - h\ij\^ij + (2.1) 



V X V X a = 47r(J„ + J, + J), (2.2) 
where the normal current J„ is given by 

J„, = ere = a(— V0 — S^a), (2.3) 
and where the supercurrent is given by 



he* 



^ (2.4) 

2mi m 

In these equations 7 is a dimensionless order parameter relaxation time; e* and m are the 
charge and mass of a Cooper pair; \a\ = ao(l — T/Tco), with Tco the mean-field transition 
temperature; a is the conductivity of the normal phase. We have also included two noise 
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terms, 9 and J, which are choosen so as to generate the correct Boltzmann weights in 
equihbrium ||T^. We then have {6) = (J) = 0, and 



{9*{x,t)9{x',t')) = 2h-fkBT6^'^\^ - ^')6{t - t'), (2.5) 

{Ji{x,t)Jj{x.',t')) = 2akBT6^'^\^-x.')6{t-t')6ij, (2.6) 

with all cross correlations zero (the brackets denote an average with respect to the noise 
distribution). In terms of the parameters in the TDGL equations, the correlation length C, = 
h/ (2m|a|)^/^, the penetration depth A = [m6/47r(e*)^|a|]^/^, the Ginzburg-Landau parameter 
K = A/^, and the thermodynamic critical field He = (47r|ap/6)^/^. 

In deriving a set of sharp interface equations it will be useful to work with a judiciously 
chosen set of dimensionless variables. Sharp interfaces between the superconducting and 
normal phases will be produced when the coherence length is small; i.e., when \a\ is 
large (this implies that we are far from the normal/superconducting phase boundary). We 
will therefore introduce a small parameter e such that \a\ = a/e^, with a a fixed constant. 
To recast the TDGL equations into dimensionless form we introduce the following primed 
dimensionless variables, 

X = Ax', t = {n-i/a)t\ = (|a|/6)i/2^' 

a = V2HM', = (a/e*7)0' , 
along with dimensionless conductivity and temperature variables 

(j = ATTK\h/2m^)a, f = kBT/[{H^/ATT)X% (2.8) 

Here we have defined A = [m6/47r(e*)^a]^/^, f = /i/(2ma)^/^, and He = (47ra^/6)^/^, so that 
^ = and A = eA. Therefore the e — limit is equivalent to taking ^, A ^ while holding 
K fixed. In these units, the TDGL equations become (we will henceforth drop the primes) 

e\dt + i<j))i) = e^{^-is^ + ^ - \^\^^ + e^9, (2.9) 
e^V X V X a = e^a (--V4> - dta) + - '^Vip*) - l^-pa + e^J, (2.10) 



(0*(x, t)0(x', t') = 2f5('^)(x - x')5(t - t'), 



(2.11) 



(J,(x,t)J,(x',t')) = aT5('^)(x-x')5(t-t')^ir (2-12) 

To further simplify, we rewrite the order parameter in terms of an amphtude and a phase, 
ijj = f exp{ix)', in terms of the gauge-invariant quantities q = a — Vx/i^ and p = (p + dtX, 
the magnetic and electric fields are 

h = Vxq, (2.13) 



e = --Vp-dtq. (2.14) 

K 

Equating the real and imaginary parts of Eq. (|2.9|) , and noting that the total current has 



zero divergence ||T5[, we arrive at the final form of the TDGL equations: 



e%f = ( -q'f)+f-f + e'C (2.15) 



e^V X V X q = (-^Vp - ^^q) - fq + e^J, (2.16) 



e^-V ■ -Vp + 9iq - /^p = -V ■ J - -C, (2.17) 

K \K J K K 

where the noise term ( has zero mean and correlations 

(C(x, t)C(x', t')) = T5(x - x')5(t - t'), (2.18) 



with the correlations of the current noise J given by Eq. ( |2.12| ) above. Note that in these 



scaled units there are three dimensionless parameters: the Ginzburg-Landau parameter k, 
which measures the coupling between the order parameter and the gauge field, the dimen- 
sionless normal state conductivity a, which determines the rate at which flux diffuses in the 
normal state, and the dimensionless temperature T, which measures the relative strength 
of thermal fluctuations. The remainder of this paper is concerned with solving Eqs. ( [^.15[ ), 
( p.l6| ), and (|2.171 ) for a moving interface. 



III. THE SHARP-INTERFACE LIMIT 



In this section we will derive a sharp interface equation from the TDGL equations, in 
the limit that e ^ while k, ct, and T are held fixed. In order to simplify the analysis 
we will assume that all quantities are translationally invariant along the direction of the 
applied magnetic field; in particular, the magnetic field is h(x, t) = /i(x, t)z, with x = 
{x,y). Extending the results to the more general three dimensional situation significantly 
complicates the calculation without offering any new physical insights. The derivation in this 
section is in the same spirit as derivations of the sharp interface equations of solidification 
from the "phase field equations" in certain distinguished limits [p!^ -[T8[] . The idea is that far 



from the interface the magnetic field in the normal phase satisfies a diffusion equation; this 
is the "outer region." Near the interface, in the "inner region," we solve the full nonlinear 
TDGL equations perturbatively in the velocity and curvature of the interface. By matching 
the solutions in an appropriate overlap region, we will see that the inner solution provides 
the boundary conditions for the outer region. This allows us to effectively "integrate out" 
the order parameter field in favor of a diffusion equation for the magnetic field with nonlinear 
boundary conditions. While this considerably reduces the complexity of the problem, the 
remaining "modified Stefan problem" is very challenging in its own right. There has been 
considerable progress in recent years in understanding the properties of this class of models 
within the context of dendritic growth; see Refs. for reviews. 

A. The outer solution 

In the outer region we assume an expansion of the form 

/(x, t; e) = /o(x, t) + e/i(x, t) + e^M^, t) + . . . , (3.1) 

q(x, t; e) = qo(x, t) + eqi(x, t) + e^q2(x, t) + . . . , (3.2) 

p(x, t; e) = po(x, t) + epi(x, t) + eV2(x, t) + . . . . (3.3) 
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Substituting these expansions into Eqs. ( p.l5|) , (|2.16|) , and ( |2.17| ), and equating terms of 
0(1), 0(e), and 0{e^), we obtain the following two sets of solutions: 

Superconducting solution. This solution corresponds to /o = 1, /i = /2 = 0, qo = qi = 0, 
q2 = J,Po = 0, pi = -C/k, p2 = -V ■ J/k. 

Normal solution. This solution corresponds to /o = /i = /2 = 0, 



V X V X qo = a ( -^Vpo - dtqo ) + J, 



(3.4) 



V- 



a ( --Vpo - dtqo ) + J 



0, 



(3.5) 



with qi, q2, pi, P2 undetermined (note that the second equation can be obtained by taking 
the divergence of the first equation). Taking the curl of Eq. ( |3.4| ) and using V ■ hg = 0, we 
see that the magnetic field in the normal phase satisfies the diffusion equation 



adtho = V^ho + V x J. 



(3.6) 



This expansion may be carried to higher order in e; order by order, the magnetic field satisfies 
the diffusion equation (without the noise term). Therefore, we find that quite generally the 
magnetic field in the normal phase satisfies 



adth = V^h + V X J, 



(3.7) 



the corrections to which are exponentially small. 



B. The inner solution 

To solve the TDGL equations in the inner region, we first afix a set of local coordinates 
(r, s) to the interface, where r measures the distance from the interface and s measures the 
arclength along the interface. In this coordinate system, the Laplacian becomes 

2 9^ d 2 2 ^ 

^ dr"^ ~^ ^ dr ~^ (^'^) ^^2 '^'^ ^ Qg^ i'^ ^) 
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where K, is the curvature of the interface [16|. Since the coordinates now evolve in time, the 



time derivatives become 

We now "stretch out" the dimension normal to the interface by introducing the scaled 
variable R = r/e. Then, keeping the lowest order terms, we have 

9. = -ii',.^ + 0(1), (3.11) 

where Vn = —r is the velocity normal to the interface. To ensure gauge invariance it is 
necessary to rescale the vector potential so that Q = q/e. The vector potential will be 
parallel to the interface, so that Q = Q{R, s,t)t, with t the unit vector tangent to the 
interface. Finally, since we want to treat the noise terms as first order perturbations about 
the equilibrium solution, we will introduce rescaled noise terms of the form 

({R, s, t) = e2C(x, t), J{R, s, t) = eJ(x, t), (3.12) 

with correlations 

{({R, s, t)C{R', s', t')) = e^f5{R - R')6{s - s')6{t - t'), (3.13) 

{Ji{R, s, t)Jj{R', s', t')) = e^af6{R - R')6{s - s')6{t - t')6ij. (3.14) 

Then in terms of these scaled variables, the order parameter amplitude F[R, s, t] e) = 
/(r, s, t; e) and the vector potential Q{R, s, t; e) = g(r, s, t;e)/e satisfy 



J_F" -Q^F + F-F' = -e (v^ + -/C ) F' - eC + 0(6^), (3.15) 



Q" _ F^Q = -e {cxvn + K,)Q' - eJt + 0{e^), (3.16) 
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where = J ■ t is the component of the current noise parallel to the interface, and where the 
primes denote derivatives with respect to R. As before, we expand F, Q, and if = V x Q = 
€ijdiQj in powers of e: 

F{R,s,t;e) = Fo{R, s,t) + eFi{R, s,t) + . . . , (3.17) 
Q{R, s, t; e) = Qo{R, s, t) + eQi(i?, s,t) + ..., (3.18) 
H(R, s, t; e) = Ho{R, s, t) + eHi{R, s,t) + .... (3.19) 



Substituting these expansions into Eqs. ( |3.15| ) and ( |3.16| ), at 0(1) we find the equilibrium 
Ginzburg-Landau equations, 

1 



^F^' -QtFo + Fo-F^ = 0, (3.20) 



Q'^ - F^Qo = 0, (3.21) 

with Hq = Qq. The solutions corresponding to an interface between the normal and super- 
conducting phases have Fq = 1, Qq = 0, and Hq = 0, for R — oo (the superconducting 
phase), and Fq = 0, Qo ~ R/V^, and Hq = 1/V2 for i? ^ oo (the normal phase). 
At 0(e), we have 

^Fi' - QlFi - 2QqQ,Fq + Fi - 3^0^^! = -(vn + \ic) - C (3.22) 
Q'l - F^Qi - 2FoFiQo = - (cjVn + JC) Q'q - Jt, (3.23) 



iJi = Q[ + JCQq. (3.24) 

The 0(e) perturbations satisfy a set of linear inhomogeneous differential equations. The 
homogeneous versions of these equations have the solution Fi = Fq, Qi = Qq, which is 
easily seen by direct substitution (this is a consequence of the translational invariance of 
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the equilibrium Ginzburg-Landau equations). This allows us to determine the value of 
Qi far from the interface without explicitly solving this set of coupled equations. To see 
this, multiply Eq. ( p.22| ) by Fq, Eq. ( |3.23| ) by Q'q, add, and integrate from — oo to R. We 
then integrate the derivative terms by parts twice, using the boundary conditions Fl{R) = 
Fi{R) = Qq{R) = for i? — oo, and using the fact that Fq and Q'q satisfy the homogeneous 
forms of Eqs. (|3.22|) and ( p.23| ). After some rearranging, we finally arrive at 

1 



HoHi + /ego 



-/C 



,V2 



R 



dx 



H,\-Q,H', + -(F[FI,-F,F^) 



-Vn r dx I (F(;)2 + a 

J — oo I 



V2 



dxQoJ^ 



R 



dxFg. 



(3.25) 



C. Asymptotic matching 

We are now in a position to match the inner and outer solutions. First, rewrite the outer 
expansion in terms of the inner variables, and take e — while holding R fixed: 



h{eR, s, t) = hQ^eR, s, t) + ehi{eR, s,t) + . . . 

dho{r,s,t) 



ho{0,s,t) + e 



R 



dr 



r=0 



O(e^). 



This expansion must match order by order onto the inner solution in the limit R 
then have the matching conditions 



(3.26) 
oo. We 



lim /io(r, s,t) = lim Ho{R,s,t), 

r— >±0 i?,^±oo 



(3.27) 



^.^aM^M)^ dH,{R,s,t) 



*±o dr 



i?->±oo OR 



(3.28) 



lim hiir, s, t) = lim 



Hi{R,s,t)-R 



dH i{R,s,t) 

"or 



(3.29) 
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Using Eqs. (|3.23|) and ( |3.25| ), we find that at the interface the magnetic field in the outer 
region is 

1 



ho{0+,s,t) + ehi{0+,s,t) = [l - e (a,s/C + r-iy^ - r^) 



(3.30) 



where a^s and F ^ are the dimensionless surface tension and kinetic coefficient for the nor- 
mal/superconducting interface, given by 



dx 



1 



V2 



Q'o 



(3.31) 



= 2a dx 



a 



(3.32) 



Some properties of the surface tension are discussed in Appendix A, and the kinetic coeffi- 
cient is discussed in Appendix B. The noise term f] is the projection of the current and order 
parameter noise onto the interface, 



f]{s,t) = -2 dxFQ{x)Cix,s,t) + 2 dxQo{x) Jl{x,s,t) 



(3.33) 

The average of the noise is easily seen to be zero; as shown in Appendix C, the noise 
correlations are 



{r]{s, t)r]{s', t')) = 2Tr~^6{s - s')6{t - t'), 



(3.34) 



so that the fluctuation-dissipation theorem is satisfied [|1J]. The boundary condition at the 
interface for the derivative of the magnetic field in the normal phase is 



dho{r, s,t) 



dr 



r=0+ 



"71^^" - Jti'T = 0+,S,t). 



(3.35) 



Finally, setting e = 1 and returning to conventional units, we have the diffusion equation 
for the magnetic field. 



dth = DV'^h + AttDV x J, 



(3.36) 



the boundary condition for the field at the interface (denoted by the subscript i), 
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hi = He 1-77^ (o"ns/C + T'-^Vn - V 
c 

and a conservation condition for the field at the interface, 



(3.37) 



h-Vh\i = D-^HcVn + ^7rJt\i- (3.38) 



In conventional units the surface tension and kinetic coefficient are 



Having now derived a set of sharp interface equations from the TDGL equations, it 
is instructive to compare our results with the heuristic set of sharp interface equations, 
Eqs. ( |l.l| )-( pr3D , which were discussed in the Introduction, and to compare our results with 
some previous work. (1) We have chosen to focus here on the two dimensional limit, in 
order to reduce the algebraic complexity of the derivation. In three dimensions we would 
reinstate the vector character of h, and the curvature /C would be replaced by the sum of 
the principal curvatures of the interface; the expressions for the surface tension and kinetic 
coefficient would be unchanged. (2) We have assumed that the material parameters of the 
superconductor are isotropic, leading to an isotropic surface tension. One way of introducing 
anisotropy is to assume that there are different effective masses along the x and y axes. The 
analysis above could then be easily modified along the lines of Ref. [0], resulting in the 
more general anisotropic boundary condition in Eq. (|1.2|) . (3) In the conservation condition, 
Eq. ( p.38| ), the critical field He appears on the right hand side multiplying f„, whereas in 
Eq. (|1.2|) the magnetic field at the interface, hj, appears. This difference is of higher order in 
e, and has therefore been neglected in our calculation. (4) The boundary condition derived 



above, Eq. ( p.37| ), depends upon the interfacial velocity, in contrast to Eq. ( |1.2| ). Such 



velocity dependent corrections also arise in derivations of the sharp interface equations of 
solidification from phase- field models [p!^ -[T8|. As discussed below, this term will have a 
natural interpretation as a local viscous damping term in the equation of motion for the 
interface. (5) We have explicitly included the effects of thermal fluctuations in deriving 
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the sharp interface equations, in contrast to Eqs. (|1.1| ) -(|1.3| ) , which neglect fluctuations 
entirely. In fact, our Eqs. ( p.36| )-( p.38| ) bear a striking resemblence to a sharp interface 
model of solidification which incorporates fluctuations, recently proposed by Karma ||20|. (6) 
Nechiporenko used the TDGL equations to study the nucleation of the superconducting 
phase, and essentially derived the inner solution in the limit of zero conductivity, and without 
noise. This results in curvature driven dynamics for the interface, without the interesting 
diffusion driven instabilities which occur when the inner solution is matched onto the outer 
solution. 



IV. THE INTERFACE EQUATION OF MOTION 

Having now reduced the TDGL equations to a sharp interface problem, we will carry the 
analysis one step further and obtain an explicit equation of motion for the interface itself 
by integrating out the magnetic fleld. This type of analysis has been used extensively in the 
study of the solidiflcation problem 



A. The boundary integral method 

We will consider a simply-connected superconducting domain which is expanding into the 
surrounding supercooled normal phase. The normal/superconducting interface is a closed 
curve C, specifled by the position vector r{s,t). As the magnetic fleld in the normal phase 
satisfles the diffusion equation, we start by introducing the Green's function G'(x, )f:|x', t') for 
the diffusion equation, 

{-D-^dt> - V'')G(x, t|x', t') = 6^''\^ - x')(5(t - t'), (4.1) 

the solution of which is (in (i-dimensions) 

G(x, t|x', t') = D{AnD\t - t'ly^/^ exp(-|x - xf /4D|t - t'\)e{t - t'). (4.2) 

Next, we (i) multiply Eq. ( |4.1| ) by h{^^t), and integrate over (x',t'); (ii) integrate the term 
involving the time derivative of the Green's function by parts, discarding the transient term; 
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(iii) use the diffusion equation for /i(x, t) to eliminate the time derivatives. The final result 



IS 



/i(x,t) 



G(x, t|x', f) V"/i(x', t') - /i(x', t') V"G(x, t|x', t 



t+ r 

dt' / d^x' 

+47r r dt' I rfVG(x,t|x',t')V' X J(x',t') 



(4.3) 



where t"*" = t + (5, with 5 is arbitrarily small, and where the subscript on the area integrals 
denotes an integral over the normal phase. Using Green's theorem, the first integral may be 
written as a contour integral along the curve C; the second integral may be integrated by 
parts with the final result 

/i(x, t) = Ho- t dt' i ds' n' ■ [G(x, t|r', t')Vh{v', t') - h{r', t') V'G(x, t\r', t')] 

J-oo Jc 



+47r r dt' i ds' Jt{r', t')G{x, t\r', t') 

J-oo Jc 

-An f dt' I V'G'(x,t|x',t') X J(x,t' 

J-oo Jn 



(4.4) 



where Hq is the external magnetic field at the boundaries of the sample (the minus sign on 
the first integral is due to the definition of the normal vector as pointing outward from the 
superconducting phase). Using the boundary conditions at the interface, Eqs. ( |3.37| ), ( ^.381 ), 
inside the integrals, we obtain 



An ^ ' ' 



+ 



An AnD J-oo Jc 
'I dt' i ds' 



dt' j^ds' G{-K,t\v' ,t')v'^ 



An J-oo Jc 

t+ 



1 - % (a„./c' + r-i< - i 



n'- V'G(x,t|r',t' 



-H^ / dt' / rfVV'G(x,t|x',t') X J(x',t' 
-00 Jn 



(4.5) 



where v'^ = Vn{s',t'), etc. This equation determines the magnetic field in the normal phase 
once the shape and velocity of the interface are specified. This equation must also hold as 
X approaches the interface itself from the normal phase. Then evaluating Eq. ( [4.5|) on the 
interface r{s,t), and using the modified Gibbs- Thomson boundary condition, Eq. ( |3.37| ), we 
obtain after some rearranging 
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47riJ J-oo 47r 



+ r7 + F, (4.6) 
where the noise term F is 

F{v,t) = H, f dt' [ d'^x'VGir,t\x',t') x J(x',t')- (4-7) 

J-00 J7V 

The correlations of F are discussed in Appendix C. 



n ■ V'G'(r,t|r',t' 



Our final interface equation of motion, Eq. ( [4 .61 ) is highly nonlinear and nonlocal. In order 
to understand the physics contained in this equation, it is helpful to dissect it and discuss 
the different terms separately. On the left hand side, the first term provides local viscous 
damping of the interface. The second term may be viewed as a nonlocal viscous damping 
term due to the eddy currents produced in the normal phase by the moving interface. 
It is this term which is responsible for the diffusive Mullins-Sekerka instability the 
propagation of the superconducting phase into the normal phase is dynamically unstable 
at long wavelengths. On the right hand side, the first term can be derived from the free 
energy difference between the two phases, and is analogous to the "undercooling" in the 
solidification problem [lIQ]. The second term arises from the surface free energy of the 
interface. The third term is a consequence of the discontinuity of h across the interface, 
which results in an effective surface current density localized at the interface. This produces 
a retarded self-interaction of the interface; in the limit of infinite diffusion constant (see 
below) this term becomes local in time and takes the form of a Biot-Savart interaction. The 
last two terms contain the thermal fluctuations of the interface. The equation of motion looks 
like a Langevin equation with nonlocal damping and colored noise. However, the retarded 
self-interaction term on the right hand side cannot be derived from an energy functional, so 
strictly speaking our equation is not of the Langevin form. As shown below, the reduction 
to a Langevin equation is complete in the inflnite diffusion constant limit. 
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B. The D ^ oo limit 



The interface equation of motion simplifies considerably in the limit that the diffusion 



constant becomes large. In this limit the second term on the left hand side of Eq. (O), 
which is responsible for the MuUins-Sekerka instability, vanishes, as does the noise term 
F. In addition, in this limit the magnetic field satisfies Laplace's equation rather than the 
diffusion equation, and the diffusion Green's functions may be replaced by a delta function 
in time multiplied by the Green's function for Laplace's equation: 

G(r,t|r',t') = --lln(i?/i?o)5(t-tO, (4-8) 

ZTT 

where R = |r(s) — r(s')|, and Rq is some long distance cutoff. The normal gradient of G is 

. V'G(r. ty, f) = '■•W-^^^g>'^<-'' ^(t - f). (4.9) 

In this limit, the interface equation of motion becomes local in time: 

r-V(s,t) = ~ - cJnslC{s,t) - - ^[a„,/C(s',t) 

+T-\^{s' ,t) - r^{s' ,t)\ y-^^ ''-^ + (4.10) 

where t = — n x z has been used. This last term on the right hand side describes the 
self-interaction of the interface, in the form of a Biot-Savart interaction of the current which 
is flowing parallel to the interface. In general this current is not constant, but rather is a 
function of the arclength s. For the moment we will ignore this complication and assume 
that the current is constant. In addition, the in our asymptotic expansion the velocity and 
curvature corrections were treated as small perturbations; this is equivalent to assuming 
that the external field Hq He- With this in mind, we see that the first term on the right 
hand of Eq. ([4.10|) is approximately {H^ — HD/Stt. With these simplifications, the equation 
of motion becomes (using f „ = n ■ dtr) 

r-n . ft. = _ .,;c - 4 / d.- I^<-- " - ^<-:- ^ ^<-'> + (4.11) 



This equation may be written in the variational form 

r-a,r = -^^ + r,, (4.12) 

where g is the metric for the interface (see Appendix D), and where the effective interface 
Hamiltonian is 

Kff H = - ^\'^' A + a^sL + ^{dsi ds'iis) ■ t(s') HR/Ro). (4.13) 
OTT lovr^ Jc Jc 

Here A is the area enclosed by the curve C and L is the perimeter of C (the necessary 
functional derivatives are carried out in Appendix D). The first term is the free energy 
difference between the two phases, the second term is the free energy of the interface, and 
the third term is the self-interaction of the interface. As mentioned above, this last term is 
a type of Biot-Savart interaction of a current sheet, with a current per unit length of H(./Att. 
The integral is therefore one-half of the self inductance of the current sheet. 

The local equation of motion, Eq. ( |4.12| ), is identical to an equation of motion for the 



interface of two dimensional ferrofiuid domains recently proposed by Langer et al. (for an 
overview, see [0). As shown by these authors, the repulsive Biot-Savart interaction tends to 
favor an extended interface, resulting in the highly branched labyrinthine patterns which are 
observed when a ferrofiuid droplet is subjected to an applied magnetic field. Labyrinthine 
structures are also observed in the intermediate state of type-I superconductors p3|, and 



the long range Biot-Savart interaction may play an important role in understanding the 



development of these patterns [22 . 



V. DISCUSSION AND SUMMARY 

In summary, we have derived an equation of motion for the normal/superconducting 
interface by starting from a fiuctuating version of the TDGL equations. In the limit of infinite 
diffusion constant (zero normal state conductivity), these equations become local and can 
be derived variationally from an energy functional. For finite diffusion constant the equation 
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of motion becomes nonlocal, and contains an additional term responsible for the Mullins- 
Sekerka instability of the moving interface. What remains to be studied is the competition 
between the self-interaction of the interface and the MuUins-Sekerka instability, and the 
implications for pattern formation in the intermediate state of type-I superconductors |2| 



After this work was completed, results similar to those obtained in Sec. Ill have been 
reported by Chapman et al. ||2^. I would like to thank Dr. Weinan E for bringing this 



reference to my attention. 
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APPENDIX A: THE SURFACE TENSION 

In this Appendix we will review some properties of the surface tension. In addition, we 
will calculate the first order corrections to surface tension near the critical value of k = 1/ -\/2, 
leading to Eq. ( |A18|) . This result has not appeared in the literature previously, and is the 
principle new result in this Appendix. 

The canonical form of the dimensionless surface tension of the normal/superconducting 
interface is 



CTri 



OO 

dx 

-oo 



(Al) 



To show that this is equivalent to Eq. ( |3.31| ), we can use the following identity: 



^0 + + F^Ql = + {Q'o? - l- (A2) 



To obtain this result, multiply the first equilibrium Ginzburg-Landau equation, Eq. ( ^.20[ ) 



by Fq, multiply the second Ginzburg-Landau equation, Eq. ( |3.21|) , by Qq, add the two 
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equations, and note that the resulting equation is an exact differential; the constant of 
integration is determined by the boundary conditions on Fq and Qq. Using this identity in 
conjunction with Eq. ( [Al| ) produces our result, Eq. ( p. 311 ). 



An explicit calculation of the surface tension requires a numerical solution of the equi- 
librium Ginzburg-Landau equations. Lacking such numerical solutions, it is still possible to 
make some qualitative observations about the surface tension. The second set of terms in 
Eq. ( p.31| ) which involve the vector potential produce a negative contribution to the surface 
tension; for sufficiently large k these terms will dominate and the surface tension will become 
negative. Detailed analytic arguments ||2^) show that for k <^ 1/v^, o^ns ~ 2^/2/3k, so that 
in conventional units a^s ~ 1.89^(if|/87r), while for k ^ l/"\/2, o- = 4(v^ — l)/3, so that in 
conventional units a^s ~ 1.10A(if^/87r). 

As first shown in the landmark paper by by Ginzburg and Landau p8|, the surface 



tension vanishes at k = 1/v^- Here we will extend these results somewhat, by expanding 
the surface tension about the k = l/\/2 limit. We will introduce a small parameter e such 
that l/(2/t^) = 1 + e. We then expand the equilibrium order parameter Fq and vector 
potential Qq in powers of e, 

Fo(x; e) = Foo(x) + eFoi(x) + . . . , (A3) 

Qo{x;e)=Qoo{x)+eQoi{x) + ..., (A4) 

where Fqi and Qoi and all of their derivatives vanish at ±00. Substituting these expansions 
into the surface tension, and noting that the 0(1) term vanishes, we obtain 

dx 2(F^o)' + 4Fooi^oi + 2g[,oQoi • (A5) 

-co '- -' 

Below we will solve the equilibrium Ginzburg-Landau equations for Fqq, and derive an 
identity for the second and third integrals above. 

The 0(1) terms satisfy the equilibrium Ginzburg-Landau equations with 1/k^ = 2: 

2F^; - Ql,Foo + Foo - F^, = 0, (A6) 
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Qm - FmQm = 0. (A7) 



Since the surface tension vanishes at 0(1), we have [p7|,pS 



-^00 — --^-^ooQoo- (A8) 

This equation may be used to ehminate Qoo in Eq. (|A6| ) , so that at k = 1/^/2 the order 
parameter amphtude satisfies 

2Ko - 2^ + Foo - F^, = 0. (A9) 

-Too 

As Fqq > 0, we can introduce a new function u{x) through Fqq{x) = exp[— 'u(x)], so that u 
is the solution to 

2u" + -1 = 0, (AlO) 

with u, n' ^ as X ^ — oo. This equation may be integrated once; using the boundary 
conditions, we obtain 



1 l\^/2 

„'=U, + V2"-^) . (All) 



2 2. 

This equation could be integrated again to obtain u{x), and therefore Fqq{x), but this will 
not be necessary for our purposes. 
At 0(e), we find 

2-^01 — Qoo-^oi — 2Q00Q01-P00 + -^01 — '^FqqFqi = — 2Fqq, (A12) 

Qoi ~ -^o^oQoi — 2F00-F01Q00 = 0. (A13) 

As in Sec. lll.B, we see that the 0(e) perturbations satisfy a set of linear, inhomogeneous 
differential equations. The method of solution is identical: we multiply Eq. ( [A12| ) by Fqq, 
Eq. (|A13D by Q'qq, add, and integrate from — 00 to x. We then integrate the derivative terms 



by parts twice, and use the fact that Fqq and Q'qq solve the homogeneous versions of Eqs. 
( |AT^) and ( [M^ ). The final result is 
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2(-Poo-^oi - -^00-^01 ) + QooQoi ~ Qoo<5oi - 



(A14) 



By substituting this result into Eq. ( [A5| ) , and performing a few integrations by parts, we 
find that the surface tension is 



(Tns = e/i (A15) 



where the integral Ii is given by 



/oo 
dxiF^of. (A16) 
-oo 

Since we have now expressed the 0(e) correction to the surface tension entirely in terms 
of Fqq, we can now use our results above to calculate the integral. By changing variables to 
u, Ii may be written as 

/•oo 



-oo 
oo 



2 / 

IX 1/2 







/•oo / I 1 \ ^/^ 

= 2 1 .„(„ + -e----) e- (A17) 

The integral is rapidly convergent, and a numerical evaluation produces Ji = 0.388. We 
therefore have our final result for the surface tension near k = 1/ \f2: 

0.388 (^-1) (A18) 



APPENDIX B: THE KINETIC COEFFICIENT 

As with the surface tension, an explicit evaluation of the kinetic coefficient defined 
in Eq. ( |3.32| ) requires a numerical solution of the equilibrium Ginzburg-Landau equations. 
In this Appendix we will make a few simple qualitative observations, and then evaluate the 
kinetic coefficient in some limiting cases. First, note that when a = k^, which in conventional 
units implies that l-nlia jm^y = 1 [see Eq. ( p^.8D ], the kinetic coefficient is proportional to the 
surface tension: 

23 



= n'^ans (a = k^). (B1) 

This is a consequence of a type of "duality" in the TDGL equations which has been previously 
noted in the context of vortex motion in superconductors |15|. However, there is generally 



no simple relation between the surface tension and the kinetic coefficient. Second, as with 
the surface tension, the terms in Eq. ( p.32 ) involving the vector potential yield a negative 



contribution to the kinetic coefficient. It is therefore possible for the kinetic coefficient to 
become negative for sufficiently large conductivity a. By balancing the two sets of terms in 
the kinetic coefficient, we see that it will become negative when ct/k ~ 1; in conventional 
units we have {2Tiha / m-y) ~ 

In order to calculate the kinetic coefficient in the large and small k limits, we can use 
the results of Ref . |^ . In the limit of small k the terms involving the gradient of the order 
parameter amplitude dominate, and we have 

f-i = {k < 1). (B2) 

In the limit of large k the terms involving the magnetic field dominate, with the result 

f -1 = -UV2- f a {k » 1). (B3) 
3 \ m7 J 

Finally, we can also calculate the kinetic coefficient exactly aX k = 1/ \/2 by using the results 
derived in Appendix A above. The result is 

= 0.388 f 1 - ) {k = 1/V2). (B4) 

At K = l/-\/2, the kinetic coefficient vanishes when 2Txfia/m'y = 1. This is to be expected, 
for at this value of the conductivity the kinetic coefficient is proportional to the surface 
tension, which also vanishes at k = 1/ 



APPENDIX C: NOISE CORRELATIONS 



In order to demonstrate that the noise term 77 defined in Eq. ( pj.33D satisfies the 
fluctuation-dissipation, we will flrst work with 
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f]iR,s,t) = -2 + 2 r dxQoJ't, (CI) 

and take the limit as i? ^ oo at the end of the calculation. The average of f] is zero. The 
two point correlation is given by 

{r]{R,s,tMR',s',t')) = 4 r rfxT rfx'F^(x)F^(x')(C(x, s,t)C(x', 

J—oo J—oo 

+4 r dxT dx'Qoix)Qoix')-^{Jtix,s,t)Jtix',s',t')). (C2) 

J — oo J —oo dXCLX 



Substituting the two point correlations from Eqs. (|3.12|) and ( p.l3| ), and integrating the 
second term in Eq. (|C2|) by parts twice, we obtain 



{f]{R, s, t)f]{R!, s\ t')) = 2T5{s - s')6{t - t') 1 2 dx (F^^feiR' - x) 

I J—oo 



+2a 



R 



dx {Q'ofe{R' - x) - Qo{R')Q'o{R')e{R - R') 



-Qo{R)Qo{RW{R' -R) + Qo{R)Qo{R')S{R - R') 



where 6 is the step function. Now take the limit i?, i?' — > oo {R ^ R'): 



(C3) 



(r/(s, t)n{s\ t')) = hm {r^{R, s, t)v{R' , s', t')) 

it. if, —> CO 



2T6{s - s')6{t - t') 1 2 r dx {F^^f 

I J—oo 



-2a lim 

R—*oo 



R 



dx {Q'or - Qo{R)Qo{R) 



(C4) 



Since lim^_ooQo(i?) = we may combine the last two terms into a single integral, and 

finally arrive at Eq. ( p.34| ), with the kinetic coefficient T^^ defined in Eq. ( p.32| ). We thus 
see that the projected noise f] satisfies the fiuctuation-dissipation theorem. 

Next, we will calculate the correlations of the noise term F defined in Eq. ( [4.7|) . Again, 
the average of F is zero, while the two point correlation is 

(F(r,t)F(r',t')) = t dh f I d^xi [ d^X2 



N 



N 



xeijekidiG{r, t|xi, ti)dkG{r', t'|x2, t2){Jj{^i, ti) J«(x2, ^2)) 
2akBTHl f dti [ rf22;^ViG'(r,t|xi,ti) ■ ViG'(r',t'|xi,ti), 

J-oo JN 



(C5) 
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where Eq. ( |2.6| ) has been used for the correlations of the current noise. By using several 
vector identities, along with the heat equation, we find that this may be written as 

{F{r,t)F{r',t')) = lLkBT[G{r,t\r' ,t') + G{r' ,t'\r,t)] 

-^2kBT f dti /rfsifii- Vi[G(r,t|ri,ti)G(r',t'|ri,^i)]- (C6) 

47r J-oo Jc 

The first term is what would be expected for colored noise, given the kernel in the second 
term on the left hand side of the equation of motion, Eq. (|4.6| ). The physics of the second 
term is unclear. It may cancel the cross-correlation (77-F), although I have been unable to 
demonstrate this convincingly. 



APPENDIX D: DIFFERENTIAL GEOMETRY OF CURVES IN TWO 

DIMENSIONS 

Here we will review a few facts regarding the differential geometry of closed curves in 
two dimensions, much of which can be found in Refs. and ||29|| . 



Let r(a;) trace out a closed curve C, with < a < 1. This parameterization is connected 
to the arclength parameterization s{a) by ds = y/gda, where g = \dr/da\'^ is the metric 
for the curve. The counterclockwise tangent vector is t(q;) = dr/da, and the unit tangent 
is t = r/y/g. In the arclength parameterization, di/ds = — /Cn, where /C is the curvature 
and fi is the outward unit normal, and dh/ds = ICt. The perimeter L of the curve C is a 
functional of r{a), and is given by 

L[r] = <j> ds 

= t da^, (Dl) 
Jo 

and the area A enclosed by C is 

^[r] = ^ (f dsr X i 
2 Jc 

- / dav{a) X T(a). (D2) 
2 JO 
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Functional derivatives of these quantities are obtained by using the Euler-Lagrange equations 
in the a-parameterization, with the results 

^ = V^/Cn, ^ = v^n. (D3) 

The Biot-Savart contribution to the energy is of the general form 

/[r] = l£dsf^ ds'i{s) ■ t(s')$(i?) 



1 da da' T{a) ■ T{a')<^{R), (D4) 

2 Jo Jo 



where R — |r(s) — r(s')|. The functional derivative of this quantity is 



where <^'{R) = d<^/dR. 
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